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Abstract 

We study invariants of three-qubit states under local unitary trans- 
formations, i.e. functions on the space of entanglement types, which 
is known to have dimension 6. We show that there is no set of six al- 
gebraically independent polynomial invariants of degree < 6, and find 
such a set with maximum degree 8. We describe an intrinsic definition 
of a canonical state on each orbit, and discuss the (non-polynomial) 
invariants associated with it. 
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1 Introduction 



The invariants of many-particle states under unitary transformations which 
act on single particles separately ("local" transformations) are of interest 
0, 1^, H, |10|, |1^] because they give the finest discrimination between different 



types of entanglement. They can be regarded as coordinates on the space 
of entanglement types (equivalently, the space of orbits of the group of local 
transformations). In this paper we study the case of pure states of three 
spin- 1 particles, or qubits. For mixed states of two qubits, it is possible 



to give a complete set of invariants ||TT|, describing the 9-dimensional space 



of orbits in terms of 18 invariants, nine of which may be taken to have 
only discrete values (for example, the signs of certain polynomials listed 
TT|). For pure three-qubit states, where the space of orbits is known 



m 



to be 6-dimensional, we can at present do no more than find a set of 
six algebraically independent invariants. We will show (Section 3) that in 
order to do this with polynomials in the state coordinates it is necessary to 
go to polynomials of order 8, and we will exhibit (Section 4) a set of six 
independent invariants; their physical meaning is discussed in Section 5. We 
will also discuss (Section 6) the possibility of finding a more convenient set 
of non-polynomial invariants. Section 2 is an introductory discussion of the 
invariants of pure n-qubit states. 



2 Pure states: general considerations 

A general theory of local invariants of mixed n-particle states has been given 
by Rains and Grassl et al. 0, |12|. Here we review the part of that theory 
that refers to pure states. 

The most general system is that of n non-identical particles A,B, . . . 
with one-particle state spaces of dimensions dA,dB,.... Let ■ i = 

1, . . . , dx} be an orthonormal basis of one-particle states of particle X; then 
the general ?7,-particle state can be written 

ijk--- 

where the sum is over values of i from 1 to d^, values of j from 1 to ds, and 
so on. By the First Fundamental Theorem of invariant theory ||15[ applied 
to U{dA),U{dB), ■ ■ ■ , any polynomial in t*-?^ " which is invariant under the 
action on |\E') of the local group U{dA) x U^ds) x ■ ■ ■ is a sum of homogeneous 
polynomials of even degree (say 2r), of the form 
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where cr, r, ... are permutations of (1,... ,r). Here tijk-. is the complex 
conjugate of t*^*^ ", and we adopt the usual summation convention on repeated 
indices, one in the upper position and one in the lower. Note that Par - is 
unchanged by simultaneous conjugation of the permutations cr, r, . . . : 

Po-T-(t) = -Po-'r'--(t) if a' — KaK"^, t' — ktk"^, 

since such a conjugation merely expresses the effect of changing the order of 
the factors in each summand in P. 

For two particles A, B there is just one permutation o", which we can 
decompose into cycles /ti, . . . , /t^ of orders /i, . . . , with /i + ■ ■ ■ + = r. 
The polynomial ^^(t) then splits into a product of polynomials P^j • ■ ■ P^si 
where P^ depends only on the order of the cycle k, which is equal to half the 
degree of P„: 



where pB = trA|^)(^| is the density matrix of particle P, with matrix 
elements 



Thus the polynomial invariants of a two-particle pure state are the sums of 
the powers of the eigenvalues of p^- These can all be expressed in terms of the 
first (Ib power-sums, which generate the algebra of invariant polynomials and 
are algebraically independent if the eigenvalues are independent. However, 
they are not independent if dA < dB, for in that case some of the eigenvalues 
of p vanish. But clearly the same argument could be used to show that the 
algebra of invariants is generated by the traces of the powers of p^, which 
is consistent because the non-zero eigenvalues of pA are the same as those 
of ps. Thus the algebra of polynomial invariants of two-particle pure states 
has a set of independent generators 



The non-zero eigenvalues of pA (or ps) are in fact the squares of the 
coefficients in the Schmidt decomposition of |^), so what we have here is the 
well-known fact that the local invariants of a pure two-particle state are the 
symmetric functions of the Schmidt coefficients. 




{Pb)1 - f'Uk. 
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3 Polynomial invariants of three-qubit states 



For the remainder of the paper we consider three spin-^ particles A,B,C. 



The classification of pure states of this system has been discussed in |1J], 
and their invariants in P, It is known ^ that the dimension of the 
space of orbits is 6; there are therefore six algebraically independent local 
invariants. We will show that there are no more than five algebraically inde- 
pendent invariants of degree less than 8, and exhibit a set of six algebraically 
independent invariants with maximum degree 8.[| 

The vector space of homogeneous invariants of degree 2r is spanned by 
functions P^rr labelled by pairs of elements of Sr, the group of permutations 
of r things. Thus there is one independent invariant of degree 2, 

Jl=Pee(t)=f%fc=(^|^) 

where e is the identity permutation, so that Si = {e}. If 5*2 = {e,cr}, the 
four linearly independent quartic invariants are 







^ ''Jljlfci'' ^1232^2 


= 


h 


= Pea{t) 


— +niifciT. . , +i2j2k2j, . , 

^ ^ii3\k2^ ^i2nki 




h 


= P<xe(t) 


iililfeiT j-i2j2k2j 

''jij2fcl'' ''i2jlk2 


= tr(p|), 


h 


= P.a{t) 


_ iiiiifciJ ^i2j2k2j 

^ii32k2^ ^i23iki 





where Pa,Pb,Pc are the one-particle density matrices: 

px = tiYz |^)(^| where {X, Y, Z} = {A, B, C} in some order. 

Thus there are at most four algebraically independent invariants of degree 
< 4. 

Higher-order invariants Pnai'^) with tt, cr G ^3 are functions of the four 
quadratic and quartic invariants if vr and a are equal or if either of them is 
the identity. To see this, note first that if vr = a, 

P (-t-") -l-hjlkl ArjrkrJ J 

J- aa\^) ■ ■ ■ ''ni<T(l)fc<T(l) ■ ■ ■ ''^rja(r)ka(r) 

= (p^)::,,(pa)-,---(pa)-, 

where r = a~^. This is a product of traces of powers of pa- But since pa is a 
2x2 matrix, the Cayley-Hamilton theorem enables us to express tr(p^) for 
r > 3 as a function of tr pA and tr p\. 

understand that similar conclusions have been reached by Markus Grassl M. 
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Secondly, if vr = e, 



= iPcfkl^,^ ■ ■ ■ {Pcfkl^^^ 

which is a product of traces of powers of pc] and similarly P-^ei^) is a product 
of traces of powers of p^. 

Thus the only sextic invariants Pt^ct which might be algebraically inde- 
pendent of the quadratic and quartic invariants are those for which tt and a 
are distinct 2-cycles, or distinct 3-cycles, or one is a 2-cycle and the other is a 
3-cycle. Moreover, in each of these categories all the possible pairs (vr, a) are 
related by simultaneous conjugation and therefore give the same invariant. 
There are therefore three possible independent sextic invariants: 

1. TT, a distinct 3-cycles, say vr = (123), a = (132). This gives 

-^5 = -P(123)(132) (t) = t^^^^^^t^'^'''^^^t^'^^'^^^thj2k^i2hkj'i3nk2 

= {PBc)fj:l{PBc)fJll{PBCyjfkl (3.1) 

where psc = is the density matrix of the two-particle system of 

particles B and C. This invariant was identified by Kempe p[ as one which 
distinguishes three-particle states which have identical density matrices for 
every subsystem. It has exactly the same form when expressed as a function 
of pab or of pAc- 

2. vr, a distinct 2-cycles, say vr = (12), cr = (23). This gives 

I5 = -P(12)(23)(t) = i*^''^'''^^*^''^''^^*^''^^'^^nj2fci^i2JiA:3^j3J3fc2 
= i.PB)]\{Pc)kl{PBc)f^k^ 

= ti[{pB® pc)pBc\- (3.2) 

3. TT a 2-cycle, say (12), and a a 3-cycle, say (123), or vice versa. These 
give 

n = P(12)(123)(t) = t'^^'^'^f^^'^'^f^^'^'^t.U2;^2^.2..fc3^.3.3fci 
= {PAC)^PA)Z{PC)1^ 

= ti[{pA® Pc)pAc\ (3.3) 

and 

/s" = ^(123)(12)(t) = f^^'^'^^f^^'^'^^t^^^'^'^tn.2;^J.2.3fcJ.3.lfc3 

= iT[{pA® Pb)pab\- (3.4) 
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Primes have been placed on the symbols for these last three invariants 
because they will not feature in our final list of independent invariants, each 
of them being expressible in terms of I5 and the quadratic and quartic in- 
variants. To show this, we write I5 in terms of 2 x 2 matrices by considering 
the 4x4 matrix pbc as a set of four 2x2 matrices Xj^ : the matrix elements 
of Xj^, labelled by (fci, /C2), are 



Then 



h = ^ixf,x^x^. 

Now we use the 2x2 matrix identity 

ii{XYZ) + ii{XZY) = tr X ii{Y Z) + tr F ii{ZX) +tiZ ti{XY) 



(3.5) 



- trXtrFtrZ 

which holds for any 2x2 matrices X, F, Z, and can be obtained by trilinearis- 



ing (or "polarising" |T5| — replace X first by X + y and then by X + y + Z) 
the cubic identity 

trX^ = ftrXtrX^ - Xf 

which in turn is obtained by taking the trace of the Cayley-Hamilton theorem. 
Apply (|3.5| ) to the matrices Xj^, Xj^^, Xj^ occurring in the expression for I^. 
The first term on the left-hand side is J5; the second is 

tr(X^X^X^) = tr(p|^) = tr(pi) 

since the non-zero eigenvalues of psc are the same as those of pA (both being 
the squares of the coefficients in a Schmidt decomposition of \^))- The first 
term on the right-hand side is 

tr(X-)tr(J^J^) = M^ipBcfS^PBcfS:, 



iPB)%{pA}Z{pAB)^ 



the second and third terms differ from the first only by permuting the indices 
ji? j2, is and therefore (after summing) are equal to it; and the last term is 



tr(X^)tr(X^)tr(X^;^)=tr(p 



B) 
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Thus ( p75|) gives 

h = 3tT[{pA (S) Pb)pab] - tr(p|!i) - tr(p|). 



(3.6) 



Similarly, using the alternative expressions for in terms of pab and pac 
gives 

/5 = 3 ti[{pB Pc)pBc] - tr(p|) - tr(p^) (3.7) 
= 3 tr[(p^ Pc)pAc] - tr(pi) - tr(p^). (3.8) 

So there are at most five independent invariants of degree 6 or less. Since 
six invariants are needed to parametrise the orbits 0, we must use at least 
one invariant of degree 8 or more. A convenient, and physically significant, 
choice is the 3-tangle identified by Coffman, Kundu and Wootters 0: 

— i-r^ — \f - f - ■ f ■ ■ f ■ ■ f, , f, , f^-^hkifi232k2fi-ink3fiijiki\^ fo q\ 

-'6 4 ' 123 \'^iii2^iiiA^JiJ2^J334,^kik3'^k2kA^ t t t | \^-'^) 

where eij is the antisymmetric tensor in two dimensions (ei2 = —£21 = 1, 
eii = ^22 = 0). The expression between the modulus signs is an SU(2)^ 
invariant (though not a U(2)'^ invariant — its phase is not invariant under 
local transformations), so its modulus is a local invariant. The invariant Jg 
can be put into our standard form of a sum of terms like ( |2.1| ) by multiplying 
the SU(2)'^ invariant by its complex conjugate 

J'JeJvhJsjeJvjSfrkskTkaksj j j j 

(where the contravariant tensor e*-' is numerically the same as e^j), and using 
the identity 

e ecd = o^d^- o^d^. 

To show that the invariants h, . . . are independent it is sufficient to 
show that their gradients are linearly independent at some point. To calculate 
these gradients in the 16- (real) dimensional space of pure states, we can treat 
f^'^ and tijk formally as independent coordinates; the fact that our invariants 
are real means that the 16 componentsof the gradient of la are the real and 
imaginary parts of the partial derivatives with respect to f^^. The results of 
calculating dla/df^^ and putting 

^000 ^ ^010 ^ ^110 ^ ^011 ^ ^100 ^ ^101 ^ ^111 ^ ^001 ^ 

T''^ = complex conjugate of f^'' 
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(where and 1 are the two possible values of i,j, k) are as follows: 

dth = (0,-1,0,1,1,1,0,1) 

dth = (-21,-81,2,8,4,10,2,8) 

dth = (0, 2 - 81, 0,6- 21, 6, 8 - 21, 2 + 21, 6 + 21) 

dth = (2 - 21, 2 - 61, 0,6- 21, 6,8- 21, 0, 8 + 21) 

= (6 - 91, 12 - 361, 6, 30 - 121, 21, 45 - 121, 9 + 61, 36 + 121) 
dth = (-8, 0, -8 + 161, 8, 8, 0, -81, 0) 

These six vectors are Indeed linearly Independent over M. 

4 Physical significance of the invariants 

The Invariant h Is just the norm of the three-party state and therefore has no 
physical significance; we will normally set it equal to 1. The three Invariants 
h,h, h are one-particle quantities, giving the eigenvalues of the one-particle 
density matrices; they are equivalent to the one-particle entropies, which 
measure how entangled each particle is with the other two together. The 
entanglement in each pair of particles and the three-way entanglement of the 
whole system are all given by the last invariant h, as follows. A good measure 
of the entanglement of two qubits A,B in a. mixed state is the 2-tangle tab 
0, which is a monotonlc function of the entanglement of formation [jl6[. The 



three-way entanglement of three qubits A, B,C in a. pure state is measured 
by the 3 -tangle 

tabc = 'Ta(bc) ~ Tab — tac 

where ta(bc) = 4detpA = 2(/i — h) is another measure (equivalent to the 
entropy of A) of how entangled A is with the pair (BC). It can be shown |^ 
that tabc is invariant under permutations of A, B and C; in fact it is equal 
to our invariant h- By solving the equations expressing the permutation 
invariance of tabc, we can now give formulae for all three 2-tangles and the 
3-tangles in terms of our invariants: 



Tab 



1 - /2 - /3 + /4 - |/, 



tac = I - h + I3 - h - \h, 
tbc = '^ + h - h - h~ \h, 

Tabc = h- 

The 3-tangle h is maximal for the GHZ state |000) + |111), whose 2-tangles 
vanish; on the other hand, h vanishes at the states p|100) -|- g|010) -|- r|001). 
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The remaining invariant, I^, is a different and independent measure of the 
entanglement of each pair of qubits. Its existence shows that the 2-tangles 
and 3-tangle are not sufficient to determine a pure 3-qubit state up to local 
equivalence. As is shown by eqs. ( |3.6|) -(|3.8| ) , this invariant is equivalent to 
any one of the two-qubit quantities kab = tr[(p^ (g) Pb)pab] (together with 
one-qubit quantities), and it relates these three 2-qubit quantities to each 
other. If we regard the hermitian operators pA and as observables, then 
kab is the expectation value of PaPb, so HAB — hh is the correlation between 
the eigenvalues of pa and pb- It is related to the relative entropy of the two- 
qubit state pab relative to the product state Pa®Pbi and is a second measure 
of the entanglement of the pair (A, 5), independent of the 2-tangle tab- 

Finally, we give the values of these invariants for some special states (all 
of which are taken to be normalised). 

For a factorised state a|lll) + 6|100), 

/2 = 1, I^ = h = a^ + h\ I^ = a^ + b^, /6 = 0. 
For a generalised GHZ state p|000) + g|lll), 

I, = I, = I,=p^ + q\ I^=p^ + q\ I, = 4p\\ 
For the minimally 3-tangled [|] state p|100) + g|010) + r|001), 

h=p' + iq' + rr, h = q'+ir' + pr, h = + {p' + q')\ 
l^=p^ + q^ + r^ + 3p'^q\^^ Jg = Q. 



5 Canonical coordinates 

An alternative type of invariant, not necessarily a polynomial in the coor- 
dinates of the state vector, is obtained by specifying a canonical point on 
each orbit. The values of the invariant functions at any point are then the 
coordinates of the canonical point on its orbit. The canonical points lie on a 
manifold corresponding to the space of orbits, and their coordinates can (at 
least locally) be expressed in terms of an appropriate number of parameters. 
One form of canonical state was suggested independently by Linden and 



Popescu 10 and by Schlienz [|T^], who pointed out that any pure state of 
three qubits can be written as 

1^) =cos^|0) (cos0|O)|O) + sin0|l)|l)) 

+ sin^|l) (r(-sin0|O)|O) + cos0|l)|l)) + s|0)|l) +^^^"11)10)) ^^'^^ 
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where 0<6',(/)<7r/4, 0<co'< 27r, and r, s, t are non-negative real numbers 
satisfying + + = 1. Simpler canonical forms, in which the number of 
non-zero coefficients is reduced to five, have since been proposed by Acin et 
al and Carteret et al 0; the latter form is 

p|100) + g|010) + r|001) + s|lll) + te'^|000) 

where p, q, r, s, t and 6 are real parameters. It is straightforward to calculate 
the invariants /i, . . . , /e in terms of either of the above sets of parameters; 
the results are not enlightening. 

We will now describe another, more intrinsically defined, form of canonical 
point whose coordinates are more simply related to Ji, . . . ,Iq. 

The three-particle state l^') has three Schmidt decompositions: 

1^) = ^ai|0i)A|$i)BC 

i 

= ^Pi\^i)B\'S>i)AC (5.2) 

i 

= ^li\Xi)c\Xi)AB 
i 

where {\0i)} and = 0, 1) are orthonormal pairs of one-particle 

states, {!$,)}, and are orthonormal pairs of two-particle states, 

the suffices indicate which of the three particles A, B, C are in which state, 
and {aj}, {(3i} and {7i} are pairs of non-negative real numbers satisfying 

al + al = al + al = 7? + ll = {^\^) = h. (5.3) 

These Schmidt coefficients, being the positive square roots of the eigenvalues 
of the one-particle density matrices pa,Pb,Pc, are related to the quartic 
invariants by 

at + "2 = tr(p^) = /2, 

+ fi\ = tr(p|) = /3, (5.4) 
it + 72 = tr(p^) = h. 

These equations have unique real non-negative solutions for a^, Pi, ■ji provided 
the invariants /i, . . . ,Ia satisfy 

/i>0, \il<h,h,h<il 

Now consider the coordinates of with respect to the canonical basis 
I'Pi) A\Gj) B\Xk)c- If the states \0j), \xk) were uniquely determined by |\E') 
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— and they almost are — then the coordinates d^'^ would be local invariants. 

However, the Schmidt decompositions do not determine the phases of l^j), 
\6j) and \xk)- We can fix these by requiring that four of the c*-'^ should be 
real: for example, we can change the phases of |(6o) and to make c^^^ and 
gioo YQoi^ then change the phases of |^o) and to make c^^^ and c"^^ real, 
simultaneously changing the phase of |xo) to keep 0°*^° and c^^^ real. (It is easy 
to show that under the six-dimensional group of phase changes of the basis 
vectors, the generic set of coordinates has two-dimensional stabiliser, so that 
the orbits are four-dimensional and therefore four phases can be removed.) 

From the Schmidt decompositions we obtain the one-particle density ma- 
trices 

PA = ^Q;-|0i)(0i|, 

i 

Ps = (5.5) 

i 

PC = '^li\Xi){Xi\- 

i 

Hence the coordinates c^^'^ satisfy 







jk 




^ ^ C C^mk 




ik 






= IkK- 







(5.6) 



To obtain a relation between the and Kempe's invariant /s, we cal- 
culate 

tr[(p^ ® Pb)pab\ 



= tr 



J 



ijk 

But 

C'''' = {(l>i\A{Oj\B{Xk\c\^) = lk{(t>i\A{Oj\B\Xk)AB- 
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Hence 



ijk 

and so, using (^.61), 

h = 3j2a^(3]\c-r-Y.^^-j:(^l (5-7) 

ijk i j 

Finally, the relation between the and the 3-tangle Jg needs a longer 
argument which we will not give here. The result is 

/6 = deti? (5.8) 

where 

i?} = {at + af)6} - Y^iPl + lfy%ki 

kl 

In order to determine how many states have the same values of the invari- 
ants /i, . . . , /e, and therefore how many further discrete- valued invariants are 
needed to specify uniquely a pure state of three qubits up to local transforma- 
tions, one would need to find the number of different sets of coordinates 
satisfying the reality conditions given above and the equations (p.6|), ( pI7| ) 
and (|5^ ), where a^, A and 7^ are determined by (|5.3| ) and (|5.4|). 
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